clear all
clc
% pressure drop line

format long;

rho=1.1397; %puknotita (37C)
v=16.7155e-6; %kinimatiki sunektikotita
D=0.0635; %diametros (m)
A=pi*D^2/4; %diatomi (m2)
L = 1; %mikos agwgou (m)
e=0.002/D; %sxetiki traxutita
Q=linspace(0,0.0267,200); %euros paroxwn
dp=zeros(200,1); %pinakas piesewn
u=zeros(200,1); %pinakas taxutitwn
coeff=zeros(200,1);

for i=1:200
    u(i) = Q(i)/A; %taxutita (m/s)
    Re = u(i)*D/v; %arithmos Reynolds
    lamda = 1/(1.14-2*log10((21.25/Re^0.9)+e))^2; % suntelestis traxutitas
    coeff(i) = (lamda*8*L*rho)/(pi^2*D^5);
    dp(i) = coeff(i)*Q(i)^2; %ptwsi piesis
end

plot(Q,dp);
axis('tight');
grid on;